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Abstract 

Three-dimensional finite element analysis was used to simulate progressive failure of a plain 
weave composite subjected to in-plane extension. The loading was parallel to one of the tow 
directions. The effects of various characteristics of the finite element model on predicted 
behaviour were examined. The predicted behavior was found to be sensitive to quadrature order, 
mesh refinement and the material degradation model. Also the sensitivity of the predictions to the 
tow waviness was studied. The predicted strength decreased considerably with increased 
waviness. More numerical studies and comparisions with experimental data are needed to 
establish reliable guidelines for accurate progressive failure prediction. 

Introduction 

Textile composites consist of interlaced tows (fiber bundles) which are then impregnated with 
a matrix material and cured. Figure 1 illustrates the architecture for a plain weave composite. The 
interlacing of the tows offers the potential for increased through-thickness strength. There is also 
the potential for reduced fabrication costs, since fairly complicated shapes can be formed using 
textile machinery. One disadvantage of textiles is the difficulty in predicting their performance. 
The complex geometry makes detailed stress analysis quite challenging. The early analyses were 
based on modified laminate theory. ( eg. References 1 ,2) In recent years there have been a few 
attempts to discretely model the fiber bundle architecture and predict internal stress states (eg. 
References 3-1 1) Reference 1 1 presented a particularly interesting progressive failure analysis of 
a plain weave composite. The results in Reference 1 1 consisted of nominal stress strain curves. 
The response of the composite was almost linear for in-plane extension and highly nonlinear for 
in-plane shear. The nonlinearity was primarily a result of progressive damage. However, little 
information was provided on damage evolution and load redistribution within the composite 
during the loading process. Also, there was no indication of the sensitivity of the predictions to 
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mesh refinement or other approximations inherent in such analyses. 

This paper has two objectives. The first is to evaluate the sensitivity of predicted progressive 
failure to quadrature order, mesh refinement, and choice of material degradation model. The 
second objective is to describe the nature of the progressive failure process for two weaves with 
very different waviness. Loading consisted of a nominally uniaxial stress along one of the fiber 
tow directions. Only mechanical loads were considered in this study. To simplify the response the 
composite was assumed to consist of an infinite number of identical unit cells in all three 
coordinate directions. 

The following sections begin with a description of the basic theory used for progressive 
damage modelling. Then the configurations will be described. Finally the results from the 
numerical simulations will be discussed. 

Theory 

There is no “right” way to model damage evolution in a textile composite that is also 
practical. It is not feasible to discretely model all of the damage, so approximation is unavoidable. 
Perhaps the simplest procedure to account for damage in a finite element model is to modify the 
constitutive matrix at the quadrature points of a numerically integrated finite element. If history 
effects are not included, the analysis of the loading becomes a series of elastic analyses. Of 
course, there are many possibilities for how to modify the constitutive matrix. Three techniques 
were used herein. The first method considered the material totally failed (ie. the entire constitutive 
matrix was reduced to essentially zero) when any allowable stress component was exceeded. This 
method will be refered to as the non-selective discount method. The second technique selectively 
reduced the rows and columns of the constitutive matrix according to the particular stress 
allowable which was exceeded. For example, if the third stress component exceeded the 
allowable, the third row was set to essentially zero to eliminate that stress component. To keep the 
constitutive matrix symmetric, the third column would also be set to zero. Zeroing the column has 
the undesirable side effect of stiffening the material with respect to the other stresses. This 
degradation method will be referred to as the selective RC method. The third technique selectively 
reduced the engineering moduli according to the particular stress allowable which was exceeded. 
Except as noted, this technique was used in the analyses. This technique will be referred to as the 
Blackketter method. 

The progressive failure analysis begins with a linear analysis of the undamaged configuration. 
Based on the calculated stresses, the initial load was scaled back so that failure would occur only 
at points which were within two percent of the maximum normalized stress. ( The stresses were 
normalized by the respective strengths.) The constitutive matrix was modified at the failure points. 
Residual forces were calculated and used to determine the incremental displacements required to 
restore equilibrium. The total displacements were updated and used to determine the new stresses. 
If no further failures occured at the current nominal strain state, the nominal strain was 
incremented to cause failure. This procedure was repeated until the nominal strain exceeded one 
percent. 

Configurations 

The fiber bundles or tows in the models were generated by translating a lenticular cross- 
section along a sinusoidal path. The waviness ratio is defined to be the ratio of the woven mat 
thickness to the wavelength. The weave consists of warp and fill tows oriented perpendicular to 
each other. In general, the warp and fill tows could be different in terms of material and shape. In 
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this study they were assumed to be identical. More details about the mesh geometry can be found 
in Reference 9. The following subsections describe the finite element meshes, the boundary 
conditions, and the material properties. 

Finite Element Meshes 

Symmetry in the material and loading was exploited so that only 1/32 of a unit cell had to be 
modeled. A wide range of mesh refinements were used, as shown in Figure 2. The crude mesh had 
only 4 elements and 42 nodes. The most refined mesh had 192 elements and 1049 nodes. The 
elements were 20-node hexahedral elements. Since only 1/32 of the unit cell was modeled, these 
refinements correspond to full cell models with 128 to 6144 elements. 

Boundary Conditions 

The periodic boundary conditions for a complete unit cell are quite simple. The appropriate 
boundary conditions for a 1/32 unit cell are a bit more complicated. Derivation of the periodic 
boundary conditions is somewhat tedious, so details will not be given here. Details can be found 
in Reference 9. The periodic conditions are listed below. Figure 2 shows the coordinate system 
assumed. 

u(a/2,y,z)= u 0 v(x,a/2,z) = constant w(x,y,c/2)= constant 

u(0,y,z) = -u(0,y,-z) v(0,y,z) = v(0,y,-z) w(0,y,z) =-w(0,y,-z) 

u(x,0,z) = u(x,0,-z) v(x,0,z) = -v(x,0,-z) w(0,y,z) = -w(0,y,-z) 


The load was controlled by specifying the magnitude of u 0 . This corresponds to uniaxial loading 
m the warp tow direction. ‘ ° 

Material Properties 

The unit cell contains two “types” of materials: the tows and the matrix pockets. Relative to 
the material coordinate system, the properties of the tows are invariant (before damage occurs). Of 
course, the properties of the tows are needed in the global coordinate system. Fourth order tensor 
transformation formulas were used to perform the required calculations. The rotation angles to be 
used in these formulas were specified at each quadrature point by using interpolation. This 
procedure was shown in References 8 and 9 to be preferable to using a single angle for the entire 
element. The particular properties used are listed below. These properties, which are 
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representative of AS4/3501-6 graphite/epoxy, are from Reference 13. 


Tow properties Matrix properties 



Modulus 

Strength 

Modulus 

Strength 

E„ 

154.27 GPa 

2342.0 MPa 

3.45 GPa 

84.85 

MPa 


10.80 GPa 

56.6 MPa 

3.45 GPa 

84.85 

MPa 

E?3 

10.80 GPa 

56.6 MPa 

3.45 GPa 

84.85 

MPa 

g 12 

7.47 GPa 

48.7 MPa 

1.28 GPa 

101.00 

MPa 

g 13 

7.47 GPa 

48.7 MPa 

1.28 GPa 

101.00 

MPa 

G 2 3 

3.33 GPa 

48.7 MPa 

1.28 GPa 

101.00 

MPa 

V 12 

0.28 


0.35 



v 13 

0.28 


0.35 



V 2 3 

0.34 


0.35 




Results and Discussions 

Most of the results in this paper illustrate the effects of characteristics of the finite element 
model on the progressive failure prediction. The effects of quadrature order, mesh refinement, and 
material degradation strategy will be considered first. Then the effect of the tow waviness on 
failure behavior be discussed. Except where indicated otherwise, the results will be presented for 
a waviness ratio of 1/3 and for the material degradation strategy in Reference 1 1, ( the 
“Blackketter method” ). 

Figure 3 shows the effect of quadrature order on the stress-strain curve. The peak stress 
obtained using 8 quadrature points (2x2x2) is 4 percent higher than that obtained using 64 points. 
The peak stresses obtained using 27 and 64 quadrature points differ by 1 percent. Damage 
initiation is predicted 3 percent earlier when 64 points are used. After the large stiffness loss 
which occurs at about 0.3 percent strain, there is even larger differences in the predictions. In 
Reference 12, non-selective discount was used for the same configuration. The difference in the 
peak stress obtained using 8 and 64 point quadratures was 10 percent. Hence, the sensitivity of the 
predictions to quadrature depends on the degradation model. The sensitivity to quadrature order 
is not suprising. For example, when, more quadrature points are used, the more extensive 
sampling is more likely to find the extremes in the stress field. One might expect a refined mesh 
to exhibit less sensitivity to quadrature order than a coarse mesh. For the meshes considered in 
this study, this was the case. Also, when failure occurs within an element, and the constitutive 
matrix is modified, the element becomes inhomogeneous. The numerical integration effectively 
fits a polynomial function to the variation of the material properties. Since the material properties 
are very different in the failed and unfailed parts of the element, it is difficult to obtain a good fit. 
In fact, there is a concern as to whether the assumed quadratic displacement functions for a 20- 
node element are sufficient to obtain a reasonable approximation regardless of the integration 
order. 

Figure 4 shows the effect of mesh refinement on the predicted stress-strain curve for the two 
waviness ratios. The 4 element model predicts the correct trends, but is quite inaccurate. The error 
in the peak stress is much worse for larger waviness ratio. The peak stress and corresponding 
strain for a coarse mesh is larger than for a refined mesh, but the peak stress and corresponding 
strain for a moderately refined mesh are not necessarily bound by the extreme cases of mesh 
refinement, (see Figure 4b) The peak stress and the corresponding strain, decrease with increased 
waviness. The ratio of the initial damage stress to the peak stress and the corresponding ratio of 
strains increase with decrease in waviness ratio. That is, not only is the strength reduced by the 
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waviness, but there is also earlier damage initiation relative to the peak stress. 

Figure 5 shows damage volume versus the nominal strain for the warp and fill tows for four 
mesh refinements. Damage volume for the resin is not shown since there is almost no damage. 
The curves are quite close for the 108 and 192 element meshes, which suggests that the curves 
might be close to convergence. As the mesh is refined, the increments in damage volume become 
smaller , but more numerous. The damage volume at the peak stress decreases with mesh 
refinement for both the warp and fill tows. This is not surprising since the corresponding stress is 
also smaller. After the peak stress is reached, there is a large increase in damage volume without 
an increase in strain (of course, the stress drops). The damage volume after this sudden increase 
is sensitive to mesh refinement for the warp tow but insensitive to mesh refinement for the fill tow. 
The curve for the 192 element mesh in Figure 5a has points labeled A and B. From A to B the 
damage initiation and growth was dominated by the inter-tow normal stress <*33. The sudden 
increase in damage at point B was due to the stress component O13 ( the xj direction is parallel to 
the fibers). 

The damage volumes just after the peak stress is reached is as follows 

Waviness Ratio Warp Fill Percent Stiffness Loss 

1/6 .51 .39 41 

1/3 .33 .35 60 

Interestingly, the damage volume is larger for the 1/6 waviness ratio, but the percentage stiffness 
loss is less. 

Figure 6 and 7 show the effect of mesh refinement and waviness ratio on damage 
accumulation during loading. The black region indicates the damage zone. The stress-strain curve 
for a particular mesh is shown above that mesh. The points labeled A,B and C indicate the 
correspondence between the strain level and the damage contours. The damage zone 
corresponding to point A indicates the initial damage. The damage zones corresponding to points 
B and C indicate the pre- and post- collapse damage states, where collapse refers to the large 
stress drop just after the peak stress is reached. Also indicated are the stress components which 
contributed to the damage. The 4-element mesh does not model the initial failure well for the two 
waviness ratios of 1/3 and 1/6. The 32 element model performs reasonably well for obtaining 
qualitative results. Further numerical studies are needed to determine how close the 192 element 
results are to convergence. For the 1/3 waviness ratio C33 dominates the initial failure. This initial 
failure appears to be an inter-tow failure resembling delamination in laminated composites. The 
collapse is characterized by 053 failure in the warp tows. For the 1/6 waviness ratio, the collapse is 
due to a significant failure of the warp tow due to 053 and cracking of the fill tow due to 022- 

Figure 8 shows the stress-strain curves obtained using the three degradation models described 
earlier: non-selective, selective RC, and the Blackketter method. The non-selective method 
predicts 22% lower peak stress than the Blackketter method. Selective reduction of rows and 
columns in the constitutive matrix results in a much larger residual stiffness after the peak stress 
than the other two methods. Selective reduction of rows and columns in the constitutive matrix 
also does not result in a large sudden drop in the stress after the peak stress is reached. 
Concomitantly, for the selective RC method there is also no sudden increase in the damage 
volume after the peak stress is reached. (This is not shown on the plot.) It should be noted that if 
variability in the strength of the constituents from one unit cell to another was included, there 
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would also be less severe jogs in the stress-strain curve. 

Concluding Remarks 

Simulation of progressive failure in a plain weave composite is extremely complex. 
Consequently, only approximate treatment is practical at this time. One of the goals of this paper 
was to examine the effect of several approximations on predicted behaviour. One obvious 
conclusion from this study is that the predictions are quite sensitive to a number of decisions 
which must be made when assembling a finite element model. Further numerical experiments and 
comparisons with experimental data are needed to establish guidelines for accurate analysis of 
progressive failure. 

Another objective of this paper was to describe the effect of tow waviness on damage 
accumulation. The results suggest that the degree of waviness not only affects the stress at which 
damage initiates, but also the type of damage which occurs. Also, the stress component 
responsible for damage changed during the progressive failure process. 
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(a) Full unit cell. 



(b) Single mat with matrix pockets removed. 


Figure 1 Schematics of plain weave composite. 
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Figure 3 Effect of quadrature order on nominal stress-strain curve. 

Waviness ratio = 1/3. Blackketter discount method was used. 
The model had 192 elements and 1049 nodes. 
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(a) Waviness ratio = 1/3 

Figure 4 Effect of mesh refinement on nominal stress-strain curve. 
Blackketter discount method was used. 
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(a) Warp tow 

Figure 5 Effect of mesh refinement on damage volume. 

Waviness ratio = 1/3. Blackketter discount method was used. 
1 indicates strain at which collapse occured. 




Nominal strain (percent) 
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Blackketter discount method was used. 







Discount method Peak Stress 
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Figure 8 Effect of discount method on nominal stress-strain curve. 

Waviness ratio = 1/3. The model had 192 elements and 1049 nodes. 
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